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This paper examines the features of a generalized position-dependent mass Hamiltonian H m in 
a supersymmetric framework in which the constraints of pseudo-Hermiticity and CIPT are natu- 
rally embedded. Different representations of the charge operator are considered that lead to new 
mass-deformed superpotentials W m (x) which are inherently ^Pi-symmetric. The qualitative spectral 
behavior of H m is studied and several interesting consequences are noted. 
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I. INTRODUCTION 



£ 
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Non-Hermitian systems admitting TT-symmetry (i.e. invariance under a combined action of parity 7 and time- 
reversal 7) have been a subject matter of intense interest [l|, 0] ■ TT-symmetry has an interesting implication that 
the whole class of Schrodinger Hamiltonians coming under its assignment namely, H = p 2 /2m + V(x) defined on 
J> , the real line x £ K, where the potential is typically V(x) = V*(—x), may possess real or conjugate pairs of energy 
eigenvalues under certain conditions related to 77 being unbroken (i.e. exact) or spontaneously broken. It has also 
been realized that the concept of TT-symmetry has its roots in the theory of pseudo-Hermitian operators and that 
pseudo-Hermiticity serves as one of the plausible necessary and sufficient conditions for the reality of the spectrum 

1- _ 

C^l . In H a set of intertwining relations 

(N ; H£ = (H ] (l.i) 



was studied in which a Hcrmitian operator £ was proposed to be expressed as a product of the charge operator C and 
parity operator 7 



C = e? (c = c + )- (1.2) 

It is straightforward to see that equations (TTTTJ) and (11.21) together imply the C77 conservation of the Hamiltonian H, 
7 being the time reversal operator 

G77H = HG77 . (1.3) 
Interestingly, it also follows from (jl.ip that the operator if it exists, also fulfills the intertwining relations 

H*(- 1 = C 1 H (1.4) 
implying that H is pseudo-Hermitian with respect to C," 1 . This can be verified as follows: 

<ip,H<p> i: -i=<i;,C 1 H<j)>=<ip,H f T 1 (f>>=< Hip,C 1 4>>=< Hip,(/)> c -i . (1.5) 
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Differential realizations for C have been considered in the literature such as C = d / d x + W (x) for the first-order [4j 
and 6 = d 2 / dx 2 + W (x) d / dx + Uq(x) for the second-order case @- The aims of such models have been to search 
for closed-form solutions and to work out the solvability criterion of the embedded Hamiltonian. 

In this article we intend to investigate these and related aspects of pseudo-Hermiticity and CTT-conservation for 
extended versions of Schrodinger equation admitting SUSY in a position-dependent (effective) mass (PDM) framework. 
The 1-D effective mass Hamiltonian H — > H m obeys (in the atomic unit defined by h 2 — 2) in real spatial coordinate 



(1.6) 



H m (x)ip n (x) = -d 



m(x) 



+ V m (x) ip n (x) = E n ip n (x), V m (x) = V m (x) + p{m) , 



where m(x) is a real valued mass function in the presence of a complex potential U m (x): 

V m (x) = V£(x) + iV£(x). 
In equation (jl.6l) . the mass-dependent function p(m) has the form Q 



p(m) 



l + bm"(x) m' 2 (x) 
— c . 

2 m 2 (x) m 3 (x) 



1 + b + a(a + b+ 1), 



(1.7) 



(1.8) 



where a and b are the usual ambiguity parameters Q typical to the effective-mass models. Position dependence in 
mass shows up in different areas of physics - semiconductors , quantum dots [l(| , 3 He clusters fll| and many more. 
A number of papers have been written on the issue of PDM in this rapidly expanding literature [12l434| . 

Note that the question of boundedness and invertibility of the operator £, assuming an explicit representation for 
it was addressed in [3] for the constant-mass case. In the PDM scenario, the problem is trickier and will be taken up 
elsewhere. 



II. PSEUDO-HERMITICITY AND efT-SYMMETRY IN A SUPERSYMMETRIC FRAMEWORK 



In the framework of supersymmetric (SUSY) quantum mechanics [35H38I ] . an underlying anticommutator K of the 
supercharges Q and Q can be explicitly constructed by specifying the following representation 



K = {Q,Q}=( CC Q * c ? c 



where Q and Q are defined in terms of the operator £ and its complex conjugate C* 

(\ — ( 



Q = 







Q 



C o 



(2.1) 



(2.2) 



Assuming polynomial expansions 

N 



N 



CC = 5> (H m ) N ~ k , C*C = £ h {H* m f- k , [J = 1 , H° m = I 2 ] 



we get by post-multiplying the first relation and pre-multiplying the second relation above by C and subtracting 

N-l N-l 



(2.3) 



\JV-fc 



k=0 k=0 

Similarly, by pre-multiplying the first relation and post-multiplying the second relation by and subtracting 



N-l 



N-l 



o = E K (H m f- k ] - E [h (H: n f- k c 



k=0 k=0 

(I2.4p and (I2.5P lead to the intertwining relations 



(2.4) 



(2.5) 



(2.6) 
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At play are also the following constraints 



Pseudo-hermiticity constraint: £^ = C =^ 6' [— x] = C[x] , 

CTT constraint: G77H m = H m G77 Q{x]H^[-x\ = H m [x\Q[x) , 

N N 

SUSY constraint: CC* = ^l k H^ k C[x]C*[-x] = J2 l k H v 



r N-k 



(2.7) 
(2.8) 

(2.9) 



k=0 



fe=0 



Finally, in the context of the 7V-th order SUSY, a mass-deformed superpotential W m (x) can be introduced which 
is given by the form [28| 



W m (x)=W(x)-y 



1 



y/m N (x) 



(2.10) 



where W (x) corresponds to the superpotential of the constant mass case. A natural consequence of (|2.10[) is that 
unlike W (a;) as in the constant-mass case W m (x) turns out to be TT-symmetric from the pseudo-Hermiticity constraint 
(|2.7[) as will be revealed below. 



III. FIRST ORDER CHARGE OPERATOR 



The first order representation of the charge operator C in a PDM scheme is given by 

e= J- j-+W(x). (3.1) 

^{x)dx V ' y ' 

From (|2.3p we have for N = 1 the projections 

tC=H m + l u eC = H; n + li- (3.2) 
Imposing the pseudo-hermiticity restriction ()2.7|) . we have the solutions: 



m(x)=m(-x), W(x) = W*(-x)-i^^-. (3.3) 
It is evident from (f2TT0l) and ((331) that 

W m (x) = TTW m (x). (3.4) 

implying W m (x) to be TT-symmetric and the mass function to be parity-invariant. As remarked earlier, W (x) 
ceases to be CPT-symmetric. A TT-symmetric W m (x) can be implemented by choosing for W (x) the form say, 
W (x) = exp(iax) + h(x) where a £ K and a non-TT h(x) can be confronted with a suitable parity-invariant mass 
function leaving W m (x) to be PT-symmetric. The following concrete example is one we have in mind 

1 7T 

W (x) = exp(iax) — sin(x), m(x) = — scc 2 (x) W m (x) = exp(iax) < x < — (3-5) 

where W m (x) is a periodic potential. 

Turning to the CTT-constraint (|2.8[) and using (|3.1|) wc get two relations. While comparison of the <9-term yields 
the difference 

~ ~ W fx) 

V m (x) - 3W m (x) - 2^^J= , (3.6) 

V m ( x ) 

the remaining part results in a first-order differential equation which can be readily integrated to provide for U m (x) 
the expression 

y/m(x) 4 m 2 16 m A 
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where A is an arbitrary constant of integration. 

A non-trivial form for V m (x) is also obtained on employing the SUSY constraint (|2.9I) for N = 1 namely 

e[x\e*[-x]=H m +h- (3.8) 

A comparison between (13.71) and (|3.8I) fixes A = — 1\. (|3.8I) is our final form of U m (x) for the N = 1 case. Note that 
the underlying CTT-invariance has the implication 

ff m e:PTV(X) = e?Ti? m v(a;) = errEip(x) = E*ey7^(x). (3.9) 

Thus if (ip, E) is an eigenpair of a CCPT-invariant PDM Hamiltonian H m , then (CTTi/;, E*) must form another eigenpair 
provided Gyit/j ^ 0. Thus CTT-invariance of ip leads to the corresponding Hamiltonian having real eigenvalues. 
From the first relation of (|2.3p and (|1.6p . the ground state in the N = 1 case has to obey 



( 4 ^W = o 



with the lowest eigenvalue — 1\. From (|3.10[) we find 



d + W*(x) 



f/>o(-x)=0, (3.10) 



exp 



y/m(y)W m (y) dy 



(3.11) 



Wo being the normalization constant. Note that ipo(x) is non-TT-symmetric. 

IV. SECOND ORDER CHARGE OPERATOR 

We now look at the following mass-dependent second order representation of the charge operator C 

' ' r W(x)^-+U (x), (4.1) 



m(x) dx 2 v J dx 
accompanied by the N — 2 SUSY representations 

CC* = Hi + hH m + l 2 I 2 . (4.2) 



as follows from (|2.3I) . For the literature on N = 2 SUSY in the constant-mass case we refer to the readers [39h44| | . 
Employing the pseudo-Hermiticity requirement (I2.7[) gives the following solutions 

m(-x) = m{x) , W m (x) = WW m (i) (4.3) 

which are similar to the N = 1 case i.e. W m (x) is TT-symmetric and the mass function m(x) is parity-invariant. Note 
that according to (|2.10[) . W m (x) is related to the constant- mass superpotential W (x) by 

W m (x) = W(x)- (-) . (4.4) 



As an illustrative example we can take this time 

W (x) = cxp(iax) — sin(x), m(x) — sec(x) a€R, < x < — (4-5) 

leading again to a periodic yT-symmetric W m (x) = exp(iax). 

Apart from (|4.3[) . the pseudo-Hermiticity condition also furnishes another relation namely 

Alto(x) = Uo(x) - TTUo(x) = W m (x), (4.6) 

which reflects the non-CPT-symmetric character of the function U.q(x) present in (|4.ip . 
Next, consideration of the CCPT requirement (|2.8I) furnishes 

— — — vn! 

AV m (x) = V m (x) - 97V m (x) = 2W' m {x) + —W m {x) (4.7) 

m 
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which is slightly different in form from the N = 1 result (|3 . T|) . In (|4.7[) V m (x) is restricted by 

V m (x) = AV m (x) + f(x) - Uo(x) + A 



where 



/(*) = 5 



mWl l (x)-—W m (x)-W' m (x) 

m. 



(4.8) 



(4.9) 



on making use of (|4.3p and (|4.6p . The constant A appears in (|4.8|) through the process of integration and is left 
arbitrary at this stage. 

In addition to (|4.7j) and (|4.8j) . the non-TT- function Ho (a;) has to satisfy the differential equation 

{•P7V m (x)}' 



m(x) 



- Uo(x)AV m (x) + W m {x) {v7V m (x)}' 



m(x) 



= 0. 



(4.10) 



2W' m (x) + —W m (x) 

771 



Substitution of (|4.8[) into (|4.10[) converts it to the form 
U' (x)W m (x)+U (x) 
which may be solved to arrive at 

U (x) 



f"(x) 



f'(x) 



+ /'(*) 



e 



W m (x) 



(4-11) 



(4.12) 
0,1,2,3,4. 

While the first two produce solutions similar to (|4.3[) . the last three respectively yields the following three relations: 

(4.13) 
(4.14) 



m(x)W m (x) m(x)W 2 n (x) mW^(x) 

being an arbitrary constant of integration. 

We now attend to the SUSY constraint (|2.9[) . Here we need to compare the five coefficients of d e 



2^-V m (x) + ^[Uo(x)+y7U (x)}-2^-W' m (x) = W 2 m (x) - (^-)'W m (x) - fc-L, 
m[x) m(x) m(x) m{x) m{x) 



1 



m(x) 



{97U (x)Y 



-^—[AUo(x)}' + W m (x)AU (x) = J—W'L(x)+[W 2 m (x)/2]', 
+ W m (x){T7Uo(x)}' + Uo(x)-?7Uo(x) = V*(x) + hV m (x) + h 



m(x) 



(4.15) 



To tackle the set of equations (|4.13l) - (|4.15l) . we observe that the second equation here can be integrated out entirely 
to have 



Alio (a;) = W m {x) +Cexp 



m(y)W m (y)dy 



where C is a constant of integration. But C has to be set equal to zero to be consistent with 
with (|4.6I) only. Incorporating it along with (|4.7I) . (I4.9[) and (I4.13|) V m (x) reads 

V m (x) = AV m (x) + f{x) - U (x) ~ h/2 , 

Then looking at (j4.8[) prompts us to identify A = — y and recast V m (x) as 

~ S 771 7T1 /i 

Vm(x) = -W' m (x) + -W ra (x) + - W 2 m (x) - U (x) - j . 



(4.16) 



So we are left 



(4.17) 



(4.18) 



We now focus on the remaining SUSY constraint (|4.15l) . This can be converted to a second-order differential 
equation 



1 hi (x) + v: 1 (-x)}' 



m(x) 



- U' (x)W m (x) + [Uq(x) - W m {x) } Uo(x) = 



K(-x) + | 



(4.19) 
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by applying the CPT-operator on both sides and rearranging. Note that the action of CPT on any function g(x) is to be 
understood in the usual sense: T7 g(x) — g*(—x) , "STg'ix) = —g* (—x) and so on. The nonlinear term Uq in (|4.19p 
is redundant and can be eliminated in the following way. Using the relation V^(—x) = 7'7V m (x) = V m (x) — AV m (x), 
(HTTSj) results in 



VZ(-x)+U (x) = f(x) 



h 



(4.20) 



Employing (|4.20l) . f|4. 19[) can be reduced to the first order form 
W m (x)U' (x) + [W' m (x) - 2f(x)]U (x) = 



(x)f'(x) 



/ 2 (z) + 



(4.21) 



Equation (|4.21j) which essentially results from the SUSY constraint (|4.15|) is consistent with the CPT equation (|4.11|) 



for Uq(x) given by (|4.12l) should we identify 9 



In terms of S = + \/l 2 — 4/2 we express Uq(x) as 



U (x) 



m(x)W 2 m (x) , W' m (x) 



^3 to (x) m"(x) 



1 



2m(x)W m (x) ' m(x) \2W m (x) 
! / r x2 



4 m 3 (x) 2m 2 (x) m(x) \2W m (a;) 
As a specific example we can go for the choice (14.5[) which would give 



(4.22) 



1 S 2 
XLo(x) = — sec(a:) exp(2iax) — — cos(a;) exp(— 2iax) 

+ — cos(x) + — sin 2 (a;) scc(a;) — - sec(a;). 
4 4 2 

Evidently Uo(x) is non-TT-symmetric. 

Let us now analyze the solution of the zero-mode equation 



i a 



exp(iax) 



(4.23) 



CV(z) = 



1 



d 2 + W*(x)d + U* (x) 



i>(-x) = o. 



(4.24) 



Two linearly independent solutions of zero-mode equation (|4.24l) may be expressed in the following compact form (see 
for details Q): 



i>j(x) = 7<jy/m(x)W m (x) exp 



where 



F 3 (x) 



m(x)W 2 m {x) + (-iyS 
2W m (x) 



F 3 {y)dy 



J = 1,2. 



(4.25) 



(4.26) 



These solutions will correspond the ground and first excited states of H mi a feature known in the quadratic SUSY 
algebra. 

Now it follows from the quadratic SUSY algebra (14.21) that the lowest eigenvalues of H m are roots of the following 
quadratic equation 



E 2 + WE + h = ^ E = - 



h+S 



h-S 



(4.27) 



It is clear that the lowest two eigenvalues 2?o,iwill be purely real if and only if the SUSY constants h,h satisfy 
following inequality 



I 2 > Al 



2 • 



(4.28) 



It may be pointed out that the condition l\ > 4?2 was identified with the reducibility of the second-order SUSY 
construction (4p| | in the context of Hermitian QM. In non-Hermitian QM, we have shown that the same condition is 
related with the reality of the spectra. 
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V. iV-TH ORDER CHARGE OPERATOR 



The charge conjugate operator G may be represented as N-th order differential operator with N coefficient functions 

jV-2 

e = ^=d N + W(x)d N - 1 + y'U j (x)d j , N = l,2,.... (5.1) 

Some of the previous results are possible to generalize. Firstly, the pseudo-Hermiticity constraint (|2.7p need to 
be compared order by order from both sides for N-th order representation (|5.1[) of C. To do this, we note that the 
contributions from the adjoint operation on the term g(x)d e may be computed using the Libneitz rule as follows 

[g(x)d^ = (-1)^ [g(x)} = (-1YJ2 [ l Vrd r {g{x)}d^] . (5.2) 

r=0 

Then order by order comparison gives the following restrictions on the coefficient functions in the charge operator 
6 given by (pTTJ) 

£ = N : 7 [m(x)] — m(x) , (mass is parity-invariant) 

i = N-l: yi[W m (x)} = W m (x) . ( superpotential is TT-invariant) 

One may compare the general result derived above with the corresponding results for N — 1 [ see (|3 . 3|) and (|3 .4|) ] and 
for N = 2 [ see (|4.3[) ]. Note that the second condition means as usual that the ReW m (x) is an even function while 
its imaginary part ImW m (i) is an odd function. 

In contrast to the superpotential W m , the functions Uj, j = to N — 2 , are not TT-symmetric. For instance, for 
£ = N-2a,nd£ = N-3we have 

U N . 2 {x)-y7[U N ^ 2 (x)] = (N-l)W' m (x), JV = 2,3,... (5.4) 



N= 1,2,3, 



(5.3) 



U N - a (x) - yi [U N _ 3 (x)} = (N —2) 
and so on. More generally, 

au n - s (x) = u N - s (x)-rm N - s (x) N c s d 



N-l I N 
2~ I ¥ 



TV = 3,4,... , 



(5.5) 



N- 



^C-iS- 1 [y7W(x)]+Y^ N - S+1 C^ [VT{U N - a+j (x)}] , 
i=i 

. (5 - 6) 

From the results, it is clear that the pseudo-Hermiticity constraints measure the amount of TT-asymmetry in the 
coefficient functions. In particular, the measure is zero for first coefficient m(x) and mass-deformed superpotential 
W m (:r). 

Next comparing the coefficients of each derivative d l for I = 0, 1, . . . N + 2 from both sides of the CCPT-constraint 
a straightforward calculation shows 



Av m (x) = v m (x) - yy v m (x) 



,N 



[2mW' m {x) + {N- l)m"W m (a;)] , N = 1, 2, 3, 



(5.7) 



Comparison for I = N — 1 gives a closed expression for the potential due to the integrability of the equation 



NV m {x) = 



C2TO 



'W m (x) + m {V^W 2 m (x) + (2N - l)W m (x) - 2Ujv- 2 (^)} 



+ 



N(N - 2) 
48 



4(2iV+ 1)1 — 1 +3iV(^-2)m( — 



+ A, N = 1,2,3,. 



(5.8) 



where we set a convention that N Cj = 0,Ujv-j = for N < j. Continuing this comparison up to the term <9°, we 
find that for all order N, only two coefficient functions in the representation of charge operator C remain independent, 
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which are the mass function m{x) and the superpotential W m . As for instance, comparing d N 2 from both sides of 
one obtains for N > 3 





i 


m 





AV m (j;) + "^Ca 



1 



U N - 2 (x) + -U' N _ 3 (x) - (N - 3) ( - ) U N _ 3 (x) 



7V+1 C 4 - 



^3 



(iV-l)(0'(-x) 



iV = 2,3,... , (5.9) 



Similar to the first and second order cases, a general nonlinear SUSY algebra can be set up. The energy in such an 
algebra are zeros of the same iV-th degree polynomial 



E N + hE 



N-l 



In-2E~ + In—iE + I 



N 







(5.10) 



from which we conclude that for an odd-order charge operator, the Hamiltonian H m possesses at least one real energy 
eigenvalue. 



VI. CONCLUSION 



In this article we have studied a generalized PDM Schrodinger equation in a non-Hermitian framework. We have 
proposed new differential realization for the charge operator and sought for the solvability of the model. Several 
interesting consequences due to PDM and non-Hcmiticity of the Hamitonian are derived. It should be noted that not 
all the results of the constant-mass non-Hermitian system are carried over to the PDM case. In constant-mass case, we 
showed that the superpotential W (x) had to be TT-symmetric to preserve CTT-symmetry and pseudo-hermiticity. In 
contrast, in the present work we have shown that the superpotential W (x) loses its TT-symmetric property. Instead 
a new mass-deformed superpotential W m (x) can be defined which turns out to be CPT-symmetric. Our work uncovers 
a new class of potentials V m {x) admitting CTT-symmetry in PDM non-Hermitian systems. We have also obtained 
extension of some of our results to a general iV-th order charge operator wherein the mass function remains even and 
the mass-deformed superpotential TT-symmetric. 
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